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Abstract 

Self-consistent mean field methods based on phenomenological Skyrme effective interactions are 
known to exhibit spurious spin and spin-isospin instabilities both at zero and finite temperatures 
when applied to homogeneous nuclear matter at the densities encountered in neutron stars and 
in supernova cores. The origin of these instabilities is revisited in the framework of the nuclear 
energy density functional theory and a simple prescription is proposed to remove them. The 
stability of several Skyrme parametrizations is reexamined. 
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I. INTRODUCTION 



The self-consistent mean-field method with Skyrme effective interactions has been very 
successful in describing the structure and the dynamics of medium- mass and heavy nuclei |l| . 
These interactions have been also widely applied to the description of extreme astrophysicai 
environments such as neutron stars and supernova cores. Actually very soon after Skyrme J2] 
introduced his eponymous effective interaction, Cameron [3] applied it to calculate the struc- 
ture of neutron stars. Assuming that neutron stars were made only of neutrons, he found 
that their maximum mass was significantly higher than the Chandrasekhar mass limit. His 
work thus brought support to the scenario of neutron star formation from the catastrophic 
gravitational collapse of massive stars in supernova explosions, as proposed much earlier by 
Baade and Zwicky |4.]. The interior of neutron stars is highly neutron rich but contains also 
a non-negligible amount of protons, leptons and possibly other particles. However micro- 
scopic calculations in uniform infinite nuclear matter using bare nucleon-nucleon potentials 
have been usually restricted to symmetric nuclear matter (SNM) and pure neutron matter 
(NeuM). Even though effective interactions are phenomenological, they can provide a con- 
venient interpolation of realistic calculations to determine the equation of state of neutron 
star cores. Mean-field calculations can be easily extended to finite temperatures and can 
thus be also used to describe the hot nuclear matter found in supernova cores and protoneu- 
tron stars. Moreover, the mean-field method allows a consistent and tractable treatment of 
both homogeneous matter and inhomogeneous matter (e.g. neutron star crusts [5]) with a 
reduced computational cost. This opens the way to a unified description of all regions of 
neutron stars and supernova cores [6|. 

Nevertheless the application of these effective forces to nuclear matter at high densities 
has been limited by the occurrence of spurious instabilities {?], In particular, Skyrme 
forces predict a spontaneous transition to a spin-polarized phase when the density exceeds 
a critical threshold which depends on the isospin asymmetry |oLJl3|] . Besides, it is found 
that for some forces the energy density of the spin-polarized phase decreases with increasing 
density. In this case, the phase transition is accompanied by a catastrophic collapse 14(, 
which is contradicted by the existence of neutron stars (note however that observations 
alone do not exclude the possibility of a ferromagnetic core inside neutron stars, see for 
instance Refs. [l5|,[l^]). Moreover, the critical density predicted within the Skyrme formalism 
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generally decreases with temperature due to an anomalous behavior of the entropy, which 
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lis instability 
M which 



is larger in the spin-ordered phase than in the unpolarized phase 
can strongly affect the neutrino propagation in hot dense nuclear matter 
is believed to play an important role in the supernova explosion mechanism and in the 
evolution of protoneutron stars 21] . However, no such spin-polarized phase transition is 



found by microscopic calculations using realistic nucleon-nucleon potentials. Indeed several 



calculations based on different methods, such as the 
method 
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26| . the Brueckner-Hartree-Fock method 



owest-order constrained variational 
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Monte Carlo method 



291 ]. the auxiliary field diffusion 
30| and the Dirac-Brueckner-Hartree-Fock method [3JJ, show that 



nuclear matter remains unpolarized well above the nuclear saturation density po both at 
zero and finite temperatures. 

The prediction of spin-ordering in nuclear matter is one of the main deficiencies of the 
mean-field method with effective forces. Different extensions of the standard Skyrme force 
have been recently proposed in order to prevent these phase transitions at zero tempera- 



ture 
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13]. In this paper, the origin of the spin and spin-isospin instabilities is revisited 



in the more general framework of the nuclear energy density functional (EDF) theory (see 



for instance Ref. 



321 ] for a review) and a simpler prescription is proposed to ensure stability 



of dense nuclear matter for any degree of spin and spin-isospin polarizations and for any 
temperature. The paper is organized as follows. The Skyrme functionals that we consider 
here are defined in Section [Til Section II I II is devoted to the discussion about the stability of 
nuclear matter. Several Skyrme functionals are reexamined in Section IIVI 

II. SKYRME FUNCTIONALS 

The nuclear EDFs that we consider here are of the form 



E = E kin + E C oui + Ei 



Sky 



where i^kin is the kinetic energy, -Ecoui is the Coulomb energy and E^ y is a functional of 
the local densities and currents (q = n,p for neutron, proton respectively): the density p q , 
the current density j q , the kinetic density r q , the spin density s q , the spin kinetic density 
T q and the spin-current tensor J qtfJiU (see for instance Ref. [l| for precise definitions). It 
is convenient to introduce the isospin index t — 0,1 for isoscalar and isovector quantities 
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respectively. Isoscalar quantities (also written without any subscript) are sums over neutrons 
and protons (e.g. Po = p = p n + P P ) while isovector quantities are differences between 
neutrons and protons (e.g. p\ = p n — p p ). The Skyrme functional E Sky is then given by 

E Sky = fd 3 r£ Sky (r), £ Sky = ]T (£f en + £? dd ) , (2a) 

J t=0,l 

£r n = Ctp\ + Ct P ptAp t + Clptn + C? J PtV ■ J t + C( Jt^Jt,,u , (2b) 

£f A = C s t s 2 t + C? s s t ■ As t + Cjs t ■ T t + Cifi + C^s t ■ V x j t . (2c) 

The spin current vector is defined by J tK = J2fny £ n^vJt^u, where e Kia/ is the Levi-Civita 
tensor. The so-called "time-even" part £° ven ("time-odd" part £° dd ) contains only even 
(odd) densities and currents with respect to time reversal. 

The coupling "constants" C£ and Cf generally depend on the isoscalar density p = p n + P P 
as follows 

C p t =a p t +b p t p a , (3a) 



C\ = a s t + b s t p a 



(3b) 



Moreover, local gauge invariance 33|, |34| imposes the following relations 



't i 



t > 



(4) 



Historically the type of functionals given by Eqs. (l2ap - (l2"c]) were obtained from the Hartree- 
Fock approximation using effective zero-range interactions of the Skyrme type l[ 

v itj = t (l + XoP^irij) + ~tj(l + a;iP CT )i [p% K r a) + <K r y) P%] 
+t 2 (l + x 2 Pa)-^Pij ■ S(rij)Pij + ^3(1 + x 3 P a )p(r) a 8{rij) 



+ -2 W (&i Pij x 5 (r y ) Pij , 



(5) 



where 



r { - rj, r 



(ji +Tj)/2, Pij = —ih(Vi — Vj)/2 is the relative momentum, P a 
is the two-body spin-exchange operator. The relations between the coupling constants in 
Eqs. ( I2bj) and (!2c|) and the parameters of the effective force in Eq. ([5]), can be found for 
instance in Appendix A of Ref . Q . 



Kutschera and Wojcik [14| pointed out that for some Skyrme forces not only is the 
ground state of NeuM polarized, but also the energy density of polarized NeuM decreases 
with increasing density. However such a catastrophic ferromagnetic collapse is ruled out by 
neutron star observations. The origin of this singular behavior can be traced back to the 
parameters t% and x-i of the Skyrme force. In particular, the authors of Ref. [ljj found that 
in order to prevent a ferromagnetic collapse of NeuM, the parameters of Skyrme forces must 
satisfy the following inequality 



t 2 {l+x 2 ) > 



(6) 



This constraint was taken into account to construct the Saclay-Lyon Skyrme parametriza- 
35] , which were fitted with the parameter x 2 — — 1 • These forces which were specifically 



tions 



developed for astrophysics, have been widely used in neutron star studies. However it has 



been found that these 
matter 
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orces predict various transitions to spin-ordered phases in nuclear 
18[ even though Eq. (J6]) was enforced. Actually this is a general 
feature of standard Skyrme forces [?], El • We will now reexamin this issue in the framework 
of the nuclear EDF theory. 



III. STABILITY OF NUCLEAR MATTER 

Let us consider the case of static uniform (possibly polarized) infinite isospin asymmetric 
nuclear matter. The Skyrme energy density, Eq. (12 ap. thus reduces to 

£s ky = Yl ( C 'Pt + C tPtn + C s t s 2 t + Cjs t ■ T t ) . (7) 

t=0,l 

Let us choose the spin-quantization axis so that the only non-vanishing components of the 
spin density s q and the spin kinetic density T q are along the z-axis. For brevity we will simply 
write s q and T q instead of s qz and T qz . In the following we will neglect the anisotropies 



induced by the polarization of matter 



Introducing the density p qa of nucleons with 



spins a =f, 4- and the kinetic density of polarized nucleons defined by 

r q . = l^f"pf , (8) 
the spin density and the spin kinetic density can now be expressed as 

Sg = p qt - Pgi , (9) 



(10) 



In fully polarized NeuM with all spins up (p = p n -\)i Eq. ([7j) reduces to 

h 2 



c-pol 

NeuM 



2M ri 



+ (C* + Cl + C& +C()p 



M + (Co + Cf + C S + CDp 1 . (11 



If the energy density is calculated from a Skyrme force in the Hartree-Fock approximation, 
we find 

C P + Cf + C s + Cl = , (12) 

CJ + C[ + C T + C 1 T = ^ 2 (l + x 2 ) , (13) 
so that Eq. (ITTj) reduces to 

T fi 2 1 1 

r„ t • (14) 



cpol 
NeuM 



& 1 /, 

+ -t 2 (l + X 2 )p 



2M„ 2 



Eq. (I12p is a consequence of the Pauli exclusion principle and the zero range of the Skyrme 
interaction. Actually as will be shown elsewhere, Eq. ( TT2l) must still be satisfied for nuclear 
EDFs that are not obtained from effective forces in order to prevent self-interactions. The 
constraint of Kutschera and Wojcik [ijj], Eq. ([6]), can thus be more generally written as 

C T Q + C{ + C T + Cf > . (15) 

If this inequality is violated, £n°uM decreases with increasing density thus leading to a 
ferromagnetic collapse. 

It is instructive to rewrite Eq. (fill) as 

P oi ft 2 

^NeuM = 2M* ^ ' 

where we have introduced the effective mass of a nucleon in a spin state a defined by 

%2 98 h2 -±[s(CZ-C?)+2s q CT} , (17) 



2M* a dr g(T 2M q 



with +(— ) for spin up (spin down respectively), and M* is the usual effective mass given by 

Ti 2 dS h 2 



-.2 ac t2 



2M* dr q 2M q 



+ [{CI - Cl)p + 2p q C\\ 
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It can be easily seen that in fully polarized NeuM, the effective mass reduces to 



2M 



2M„ 



2M 



2M r 



- + (C+ +C{)p 
+ t 2 {l + x 2 )p 



(19) 



so that Eq. ( 1161) coincides with Eq. (JHJ). Setting x 2 = — 1 as in the Saclay-Lyon Skyrme 
forces 35( therefore implies that the effective mass of polarized neutrons is equal to the bare 
mass. 

We have seen that the constraint of Ref. 14j is equivalent to the requirement that the 
effective mass of polarized neutrons remains always positive. However the ground state of 
NeuM (and more generally that of isospin asymmetric nuclear matter) could still be polarized 
as shown below. 



A. Landau stability criterion 

The stability of unpolarized homogeneous nuclear matter with respect to spin and spin- 
isospin polarizations has been generally addressed using the Landau Fermi-liquid theory (see 
e.g. Ref. {37]]). In this theory, the elementary excitations of the liquid at low temperatures 
are described in terms of quasiparticles which are in one-to-one correspondence with single- 
particle states of the non-interacting Fermi gas. Any small change 8n(k) in the distribution 
function of quasiparticles with wave vector k leads to a change 88 in the energy density, 
which can be expressed (up to second order) as 

5£ = Swf e{k)5m + llwflwr v{k ' k '^ m w (20) 

where e(k) is the energy of a quasiparticle with wave vector k and v(k,k') is the residual 
interaction between quasiparticles with wave vectors k and k'. 

In pure NeuM, the residual interaction (neglecting tensor interaction) can be expressed 



as 

v NcnM (k,k') = 1 [F NeuM (k,k') + G NeuM (k,k')a-a'] (21) 
where N is the density of states at the Fermi surface given by 

, (22) 

n it 2 
7 



with k F = (37r 2 p) 1 / 3 . We have also introduced the Pauli matrices & and a' in order to take 
into account the spin of the quasiparticles. Small perturbations involve only quasiparticles 
at the Fermi surface, i.e. with k = k' = kp. We can thus expand each term in the residual 
interaction in Legendre polynomials Pe(cos8) where 9 is the angle between k and k 1 . For 
instance, 

+00 

F NeuM {k, k') = euM ^(cos 6) (23) 

where i^ NeuM are dimensionless Landau parameters. Similarly, we can define Landau param- 
eters G^ euM . For the Skyrme functional, the only non-zero Landau parameters are of order 
I = and £ — 1. The stability of the initial state is ensured if any change in the energy per 
particle e = S/p is positive. This condition leads to Landau's criterion 



^NeuM 



>-{2£ + l) , (24a) 



GT uM > -(21+1) . (24b) 

In particular, the condition on G^ enM guarantees that NeuM is stable against small fluc- 
tuations of the (isoscalar) spin polarization I a = s / p = (p-j- — p^). This can be seen by 
expanding the energy per particle up to second order in I a 



e{I a ) ~ e(0) + 



ld 2 e 



2dl, 
with 

d 2 e 



K (25) 

/<T=0 



on 



^(l + G^ cuM ) . (26) 



3M, 

The first order term vanishes because of the requirement that the unpolarized phase be an 
equilibrium state. 

Using the Skyrme functional, we find 



G NeuM = 2N 



c° + c s 1 + k 2 F (c?; + cf) 



Now if the Skyrme functional is fitted to a realistic equation of state of NeuM 
that Cq + C{ < 0, which according to Eq. ( fl2i) implies that 



(27) 
we find 



Co + G{ > . (28) 



Ferromagnetic instabilities are therefore mainly due to the coupling constants Cj . In order 
to fulfill the Landau's stability condition C^ euM > — 1 at any density, we must have 1 

C T + Cf>0 . (29) 

The absence of a ferromagnetic transition in NeuM does not generally forbid the occurence 
of spin-ordered phases in asymmetric nuclear matter. Let us consider in particular SNM. 
The most general form of the residual interaction (neglecting tensor interaction) can be 
expressed as 

v s ™(k,k') = [F(k. k') + F'(k. k')T ■ t' -i- G(k. k')a ■ a' + G'(k. k')a ■ a'f ■ t'] (30) 
where N is the density of states at the Fermi surface given by 

N^'-f^ , (31) 

ft 7l Z 

with kp = (37r 2 p/2) 1//3 and M* is the isoscalar effective mass defined by 

M _ 2M 2 1 1 

~M* ~ ' M~W n + W p ■ {62} 

We have also introduced the Pauli matrices f, f in order to take into account the isospin 
of the quasiparticles. As before, we can define dimensionless Landau parameters Fi, F£, Gi 
and G' e . The Landau's stability conditions are in this case given by 

F e >-(2£+l) , (33a) 

F' t >-{2t+\) , (33b) 

G e >-(2£+l) , (33c) 

G' e >-{2£+l) . (33d) 
The Landau parameters F and Fq are related to the usual compression modulus 

K v = ^{l + F ) , (34) 



1 This inequality is not strictly required if the coefficients C t s are allowed to depend on the density according 
to Eqs. (|3a[) and (|3b[) and the term in (6g + 6f )p a dominates at high density. However for modern Skyrme 
parametrizations such a situation does not arise because 3a < 2. 
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and symmetry energy 



J 



h 2 h 2 

w (1 + F »» 



(35) 



respectively. The conditions on Go and G' ensure that SNM is stable against small fluctu- 
ations of isoscalar and isovector spin densities respectively. These Landau parameters can 
be expressed in terms of the spin asymmetry coefficient, defined by 

ld 2 e 

la 



h 2 k 2 

a ^FO . 

6M* 



" 2di 2 

and the spin-isospin asymmetry coefficient, defined by 

1 d 2 e 



1 + Go) 



(36) 



2dn 



/ar=0 



6M* 



(1 + G' 



(37) 



where I CT = S\j p = (p n ^ — p n ^ — p p ^ + p P i)/p- Using the Skyrme functional, the Landau 
parameters Gq and G' are given by 



G = 2N 



^0 + ^O^FO 



G' = 2N 



(39) 



The stability of SNM at any density thus requires 



Cf>0 



(40) 



These two conditions entail Eq. (1291) . Note that Landau's stability conditions allow one of 
the coefficients Cf to be negative provided their sum remains positive. 

Landau's stability conditions, Eqs. (I24bp . (l33cp and (133dl) . guarantee that the unpolarized 
state is locally stable (metastable) against small fluctuations of the spin and spin-isospin 
polarizations. But this criterion does not necessarily imply that the unpolarized state is the 
ground state, i.e. the state with the lowest energy. In particular, the ground state could still 
be polarized with finite values of I a and I aT . Moreover we have only considered so far the 
two limiting cases of SNM and NeuM. However, the outer core of neutron stars is formed 
of isospin asymmetric nuclear matter whose composition varies with depth. We thus need a 
more general stability criterion. 
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B. General stability criterion 



Asymmetric nuclear matter is stable with respect to any degree of spin and spin-isospin 
polarizations whenever the energy density £ po1 of the polarized state is larger than the energy 
density £ un P o1 of the unpolarized state (for a given density p). Using Eqs. ([7]), (TT7|) and (fTSl) 
we find 

gP ° l = E IrV + C V + C°(s n - s p f + C V + G{{ Pn - Pp f (41) 



</.CT 



qa 



which for unpolarized matter (i.e. s q = 0, T q = 0) yields 



2M. 

with 



^ UnP01 = E ^L T <i + °oP 2 + C i(Pn ~ Pv? , (42) 



3 



r 9 = ^(37r 2 ) 2/3 pf . (43) 



The difference can thus be expressed as 

£ P oi _ £ Unp oi = J- ^K o1 - r q ) + C s s 2 + C?(s n - s p f 

q 1 

+C£sT + C^(s n - s p )(T n - T p ) (44) 

where t^° 1 = r q ^ + r q i is the nucleon kinetic density in the polarized phase. The absolute 
stability of the unpolarized phase can be insured by requiring each term be separately 
positive so that £ po1 > £ un P o1 . Now the first term in Eq. (|44|) is always positive since 
mechanical stability requires M* > and the Pauli exclusion principle implies that t^° 1 > r q . 
Let us also remark that (s n — s p )(T n — T p ) > because r qa increases monotonically with p qa . 
The following constraints 

C s t > , (45a) 

and 

C t T >0 , (45b) 

therefore guarantee the absence of any spin-ordered phase transitions in asymmetric nuclear 
matter. It is readily seen from Eqs. (1271) . fl38|) and fl39|) that these inequalities enforce 

11 



Landau stability conditions, Eq. ( 124bj) in NeuM and Eqs. ( I33cp and (133dj) in SNM. Since 
Eqs. (145 ap and (I45bj) ensure the stability of asymmetric nuclear matter, they obviously 
prevent a ferromagnetic collapse of NeuM as can be seen from Eq. (Ti"5"j) remembering that 
Cq + C\ > as a consequence of M* > 0. 



C. Anomalous behavior of the entropy 



We have seen that the stability of nuclear matter requires that Cj > 0. However, these 
coefficients cannot take arbitrary values. From Eq. (14"]) , large positive values of Cj translate 
into large negative values of C/ which, in certain circumstances, can lead to instabilities in 
finite nuclei whose consequence is a major rearrangement of the single-particle spectrum 441 ]. 
We will now show that these coupling constants can be further constrained by requiring the 
stability of nuclear matter with respect to any degree of spin and spin-isospin polarizations 
at non-zero temperatures. 



It was shown in Refs. 



17 



181 ] that not only do Skyrme forces predict a ferromagnetic 



transition in NeuM above a certain critical density, but worse this density decreases with 
increasing temperature due to an anomalous behavior of the entropy. This argument can be 
easily transposed to asymmetric nuclear matter as follows. At low temperatures (compared 
to nucleon Fermi energies), the difference between the entropy density iS po1 of the polarized 
state and the entropy density ^"p 01 of the unpolarized state is approximately given by 



^pol ^unpol y 



71 



q.a 



l TM* qPq 

2h 2 kl q 



M* 

qa 



2p. 



qa 



Pq 



1/3 



(46) 



Now because the polarized phase is more ordered than the unpolarized phase, its entropy 
according to Boltzmann's definition should thus be lower, i.e. iS po1 < 5 unpo1 as found in 



realistic calculations 
find from Eq. (T4"6]) 



23 



-25. 



29|. Since this should be true for any isospin asymmetry, we 



This condition reduces to that of Ref. 



M qa fPqc\ < 2 2/3 



Pq 



(47) 



17J in the limiting case of fully polarized NeuM. 



Equation ( 14"7|) can be equivalently expressed as (q' ^ q) 



EL 



I V/ 3 
- TI 

aq J- ± i 



+ 



aq' 



I -EI, 



I ) 1/3 

1 aq) 

TI, 



< 2 



(48) 



aq 



aq' 
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With I aq = (ptf - Pql)/p q , 

2M* 

E = (C^ + C^)p ql ^ , (49) 
2M* 

r = {cT-c T 1 ) Pql -^- . (so) 

We have found numerically that the inequalities (j4"5j) can be satisfied for any degree of spin 
and spin-isospin polarizations, i.e. < \I a q\, \I<rq'\ < 1 5 provided 

E 1 <E<E 2 , (51a) 
T = , (51b) 



with Hi ~ —0.21 and S 2 ~ 0.54. We have also found solutions of pSjl for |T| > T C (S) > 0. 
But it can be seen from Eq. ( 150]) that such solutions cannot exist for all densities and must 
therefore be excluded. Inserting Eq. (j4"9~j) in Eq. ( 15 lap using Eq. f flBj) yields 

P A^ + Cl)-E 2 {Cl + Cl)]-E 2Pq ,{Cl-Cl)<E 2 ^- , (52a) 

p q [(C^ + C^-E 1 (C^ + Cl)}-E lPq/ (C T -Cl)>E 1 ^ I - . (52b) 

The terms in p q > always satisfy the above inequalities. This is a consequence of the positivity 
of M* for any density and isospin asymmetry which requires that Cq+C{ > and Cq — C\ > 
0, as can be seen from Eq. (fT8j) . The conditions (I52al) and (j52bp can be ensured for any 
density p q by imposing that the associated terms be respectively negative and positive 
leading to 

^i{Cl + Cl)<CT + CT<~ 2 {Cl + Cl) . (53) 
On the other hand, Eq. (I51bj) implies 

= . (54) 
Combining these inequalities with Eqs. (j45b]), we arrive at the following restrictions 



CZ = Cf, 0<C? < l ~Z 2 {CZ + Cl) . (55) 

Eqs. fl55|) guarantee that asymmetric nuclear matter remains unpolarized at finite tempera- 
ture T since the free energy density of the polarized phase, defined by J- 1301 = £ po1 — TS po1 , 
is always higher than the free energy density J^ un P° l = £ un P o1 _ j^unpoi Q f ^ e un polarized 
phase. 
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IV. STABILITY OF SKYRME FORCES REVISITED 



Conventional Skyrme forces have been shown to predict various spin and spin-isospin 
instabilities in nuclear matter (7, 9, 10, , |l8]. We have seen in the previous section 

that for a nuclear functional given by Eqs. (l2aj) ; (12bp and ( 12c)) the stability of asymmetric 
nuclear matter at any temperature can be ensured by imposing Eqs. (145 al) and (I55I) [the 

14|, Eq. ([6]) and more generally Eq. ffT5]) . prevents a collapse 



constraint proposed in Ref. 



of polarized NeuM, but does not forbid a ferromagnetic transition]. While the coefficients 
Cf are generally positive (at least for not too high densities), standard Skyrme forces yield 
negative values of at least one of the couplings constants Cj . The origin of the instabilities 
can therefore be traced back to the time-odd terms St ■ Tt, which are related to the time- 
even terms Y\ v Jt^uJt^u due to gauge invariance @. Since the seminal work of Vautherin 
and Brink [38], it is commonly taken for granted that the spin-current tensor (which is 
usually approximated by the spin-current vector J q ) is small in nuclei, and most Skyrme 
parametrizations therefore neglect them. We have tested this assumption by computing 
the HFB energies with and without the J 2 and J 2 terms (denoted respectively by -Ehfb 
and -Ehfb) f° r an even-even nuclei with Z,N > 8 and Z < 110 lying between the proton 
and neutron drip lines (Note that when the J 2 and J 2 terms are included, the associated 
time-odd terms in Cj play a role in the exact treatment of the masses of odd nuclei, but not 



in the equal- filling approximation [39], which we adopt here, as in all our previous papers). 
The differences AM = -Ehfb — -^hfb are shown in Fig. [1] for the Skyrme parametrization 



BSkl7 



40 



41] which was originally fitted with the J 2 and J 2 terms, and for SkI2 42) which 



was not. The impact of the J 2 and J 2 terms is quite large, reaching about 20 MeV for the 
heaviest nuclei. The impact of dropping or including the J 2 and J 2 terms is logically found 
to be correlated to the amplitude of the Cf = —C( coupling constants, especially Cq . For 
instance, in the case of the SLy4 35] interaction ( CT = —17.21 MeV fm 5 ), the HFB energy 
is affected by no more than 5 MeV, while for SkO 46] (Cq = —220.54 MeV fm 5 ) values up 
to 30 MeV can be reached. Adding or removing the J 2 and J 2 terms a posteriori without 



refitting all the parameters of the force can thus lead to significant errors. However in all 

previous studies of spin and spin-isospin instabilities in nuclear matter 

time-odd terms St ■ Tt were taken into account whereas the Skyrme forces were generally 



fitted without the J 2 and J 2 terms. This treatment not only violates gauge symmetry but 
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FIG. 1: Differences between the HFB energies estimated with and without the J 2 -terms for two 
Skyrme forces SkI2 (upper panel) and BSkl7 (lower panel) for all even-even nuclei with Z,N > 8 
and Z < 110 lying between the proton and neutron driplines. 



also introduces inconsistencies in the residual interaction hence in the Landau parameters 



(see the discussion in Section III of Ref. 



43| and also in Section 5D of Ref. 



We have therefore reexamined the stability of several Skyrme parametrizations for which 

45), SLy4 [35], SkIl-SkI5 3, 



the J 2 and J 2 terms were not included in the fit: SGII 
SkO 



46| and LNS [47j. The parametrization SGII 45J was constructed in order to improve 



the Landau parameters Go and G' and the description of Gamow- Teller resonances in nuclei. 



The Skyrme Saclay-Lyon forces and especially the parametrization SLy4 



351 ] . have been 



widely used not only in nuclear physics, but also in neutron-star studies because these forces 



were constrained to reproduce a realistic neutron-matter equation of state. The Ski [42] 
forces were all constrained (except for Skll) to reproduce the isotopic shifts of the root 
mean square charge radii of neutron rich Pb and Ca nuclei. Forces SkI3 and SkI4 were 
constructed with non-standard spin-orbit couplings. For the parametrization SkI5, the 16 O 
ground-state data were excluded from the fit. We have also included the parametrization 
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SkO |46| from the same group. The parametrization LNS [47| was fitted to Brueckner 
calculations. The Landau parameters in SNM and in NeuM calculated at saturation density 
Po, with and without the terms in Cj are shown in Table [H For comparison we have 
also indicated the predictions from Brueckner-Hartree-Fock calculations in SNM Q] and 
from realistic calculations based on the renormalization group approach in NeuM [49]. As 
can be seen in Table setting Cf = in Eqs. ([27}, fl33cj) and §M) tends to reduce 
the discrepancies between the different Skyrme functionals and generally leads to a better 
agreement with realistic calculations, especially for G' . In particular, the new values of G' 



lie closely inside the empirical range of 1.0 ± 0.1 deduced in Ref. [50] from the analysis of 
Gamow- Teller resonances and magnetic-dipole modes in finite nuclei. The improvement is 
quite spectacular for the parametrization Skll. In the case of LNS, setting Cq = actually 
deteriorates the value of the Landau parameter Go since the latter was directly fitted to 
the value obtained from realistic calculations. Table [TT] shows the critical densities of the 
spin-ordered phase transitions according to Landau's stability criterion. It can be seen that 
dropping the terms St ■ Tt eliminates the instabilities in almost all Skyrme forces. This 
prescription is also consistent with Eqs. ( 155]) and therefore prevents an anomalous behavior 
of the entropy thus ensuring the stability of nuclear matter for any temperatures. 

Moreover, setting Cf = is the only prescription which guarantees the Landau stability 
conditions of Eqs. (I24bl) . f!33cj) and f!33dj) at any density both for £ = and £ — 1. Indeed 
the Landau parameters G±, G[ in SNM and Gf euM in NeuM, are given by 

G 1 = -2N C%4 , (56) 
G[ = -2N Cfk 2 F0 , (57) 

Gf uM = -2NkUC£ + Cl) . (58) 

Requiring G x > -3, G\ > -3 and Gf euM > -3 for any density thus leads to Cf < 0. 
Combining these inequalities with Eqs. PU|) yields Cj = 0. Adopting these particular values 



tends to be supported by the following basic sum rules of Landau Fermi liquid theory 



5l| 



F e /(2£+l) l + F;/(2£+l) 

+ l + G,/(2£+l) + l + G'J{2£+\) ' l ° Jaj 
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and 



So 



n 



F e /(2£ + l) l + F;/(2£ + l) 

Ct , n G'l 



+ 9 : 



l + G t /(2i+l) l + G' e /{2£+l) 







(59b) 



Even though Skyrme forces generally violate these sum rules, the prescription C/ = sig- 
nificantly improves the second sum rule as can be seen in Table IHIl It is quite remarkable 
that dropping the terms St ■ T± not only removes all kinds of instabilities in nuclear matter 
but also improves the internal consistency of the nuclear functional. Nevertheless with this 
prescription, the Landau parameters Gi, G[ and G?^ euM all vanish leading to unrealistic 
effective masses in polarized matter. Indeed, according to Eqs. (jTTl) M*^ = M*^ = M* which 
obviously holds in the limit of vanishing spin polarizations but is otherwise contradicted 



by realistic calculations 23M25I. |29|, [3JJ- In particular, these calculations indicate that in 
polarized NeuM M n ^ > M n ^ whenever p^ > p±. Imposing the less stringent stability con- 
ditions (1551) leads to a splitting of effective masses but with a wrong sign. This deficiency 
calls for further extensions of existing Skyrme functionals. 

In the discussion above, we have implicitly adopted the point of view of the nuclear EDF 
theory 32j that the different terms appearing in Eqs (12bj) and (j2cl) can be a priori considered 
as independent from each other (apart from the requirements of gauge invariance). It is 
therefore perfectly legitimate to set C/ = — Cf = 0. However in the framework of effective 
forces, the coupling constants are uniquely determined by the parameters of the force. In 
particular, the coefficients and Cj are now given by 



Xq 



—t (- 

24 3 V2 



^3 p 



(60a) 



C{ 



48 



hp c 



- Xi 



+ x 2 



(60b) 
(60c) 



CI 



T 



We have therefore studied the stability of the 



(60d) 



with the J 2 and J 2 terms: SkP 



ew Skyrme parametrizations which were fitted 



52j, SLy5 [35|, SkO' |46j, SkX |53| and BSkl7 j40|,|4l|. The 
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parametrization SkP, which was specifically designed to be used both in the particle-hole 
channel and in the particle-particle channel, is still used nowadays. The forces SLy5 and 
SkO' were fitted following the same protocol as SLy4 and SkO respectively, but they include 
the contribution of the J 2 and J 2 terms. The force SkX [53] was constructed in an attempt 
to improve the description of single-particle energies. BSkl7 is the force underlying our 



nuclear mass model HFB-17, based on the Hartree-Fock-Bogoliubov method 40|,|4l|. With 



this model we were able to fit with an rms deviation of 0.581 MeV the 2149 measured masses 
of nuclei with N and Z > 8 given in the 2003 Atomic Mass Evaluation 54j, while at the 
same time constraining the underlying Skyrme force to fit properties of SNM and NeuM, as 
determined by many-body calculations using realistic potentials. The values of the Landau 
parameters in SNM and in NeuM are shown in Table IIVI and the critical densities for the 
onset of instabilities are shown in Table IVl For those few Skyrme forces which include the J 2 
and J 2 terms, nuclear matter is therefore unstable because of the tight correlations between 
the different coupling constants in the energy density. In order to illustrate the impact of 
the J 2 and J 2 terms and their time-odd counterparts on the stability of nuclear matter, we 
have plotted in Fig. [2] the difference between the energy per particle in fully polarized NeuM 
and in unpolarized NeuM for the parametrizations SLy4 and BSkl7. Both have been fitted 
to a realistic equation of state of NeuM, but BSkl7 includes the J 2 and J 2 terms while SLy4 
does not. Removing all instabilities requires that we impose Cf > and Cj = 0. Since the 
first term in £ of the Skyrme force is generally associated with the long-range attractive part 
of the nucleon-nucleon interaction while the density-dependent term in is related to the 
strongly repulsive short-range part, the coupling constant Cg can be made positive for any 
density by choosing x < 1/2 and x 3 > 0. With t < and t 3 > 0, the coefficient C{ will be 
positive, at least for not too high densities. Spin- and spin-isospin instabilities thus generally 
arise mainly from the coupling constants Cj and Cf , which in turn are generated by the 
momentum-dependent terms in t\ and ti- Using Eqs. f!60cj) and (160dl) . the conditions Cf = 
entail t\ = t 2 and x\ = —x 2 - Imposing these constraints would leave no degree of freedom 
for adjusting surface properties of nuclei, which also depend on the momentum dependent 
t\ and t 2 terms through the coupling constants C^ p . This would also have an impact on 
the coupling constants CJ which determine the nucleon effective masses, Eq. ([18]) . There is 
little doubt that such a force would yield poor results when applied to nuclei. Since thermal 
effects on the spin polarization are rather small for temperatures found in protoneutron stars 
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and supernova cores 17] . one may be tempted to require the stability of cold nuclear matter 
only. But even in this case, it was shown in Refs. that it is not possible to avoid 

spurious transitions to spin-ordered phases in nuclear matter above 2-3 times saturation 
density, and at the same time giving reasonable properties of SNM. We have found that 
the critical densities above which instabilities occur are even lower when more nuclear data 
are included in the fit of the effective interaction. In particular, conventional Skyrme forces 
fitted to essentially all experimental nuclear mass datapredict a ferromagnetic transition in 
NeuM at a density slightly above saturation density [13J (see also Table N} . 

The stability of cold nuclear matter can only be restored by including additional com- 
ponents in the Skyrme interaction, thereby inducing new terms in the energ y d ensity. Two 
different extensions have been recently proposed. Margueron and Sagawa 12| considered 



extended Skyrme forces with two new £3 like terms depending on the nucleon spin densities 
s q of the form 

|*5(1 + x s 3 P a )s(rf 6( rij ) + \t s i (l + atfP^r) 2 ^) . (61) 

In the energy density, Eqs. (12bl) and (l2c|) . these new terms modify the coefficients C t s . 
The additional parameters were adjusted so as to ensure the Landau stability conditions 
G > -1, G' > -1 and G^ euM > -1. The nuclear mass model HFB-17 0, Q was thus 



refitted with these new terms 55[. With this extended Skyrme force called BSkl7st, it 



was possible to maintain the quality of the HFB-17 mass model, and at the same time the 
Landau parameters were adjusted so as to remove the spin and spin-isospin instabilities 
present in the original force BSkl7. Unfortunately instabilities were still found for finite 
spin and spin-isospin polarizations 55[. The reason is that terms of the form given by 
Eq. ( 16TT) . do not change the coefficients Cf and consequently, Eq. ( 1441) is not guaranteed to 
remain positive for any spin and spin-isospin polarizations. Moreover, as noted in Ref. 55] 
the contributions of Eq. (IfJTj) to the energy density cancel in fully polarized NeuM so that 
BSkl7st still predicts a ferromagnetic collapse of NeuM as BSkl7 does. The extension of 
Ref. 12| does not affect the coefficients Cf hence also the effective masses of spin-up and 



spin-down nucleons are not affected, as can be seen from Eq. ( JT71) . This means that if the 
original Skyrme force violates the constraint ( H71) . this will still be the case for the extended 
version of this force. 

Alternatively, instabilities can be avoided by introducing into the force, density-dependent 
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FIG. 2: Difference between the energy per particle in fully polarized neutron matter and in un- 
polarized neutron matter for two Skyrme forces SLy4 and BSkl7, with (dashed line) and without 
(solid line) the J 2 and J 2 terms and their time-odd part. The black dots indicate the densities at 
which the difference vanishes. 



generalizations of the usual t\ and t 2 terms of the form 



13| 



+ x 4 P a )^ [p% p(rfS( rij ) + <5(r 4i )p(r)*p?.] + f 5 (l + xsP*)^Pij • P (ry5(r tJ ) PlJ .(62) 

These new terms modify the coefficients C 4 T , C 4 T , Cf p and C t As thus providing more flexibility 
to remove instabilities without deteriorating the fit to nuclear data. We have constructed a 
new nuclear mass model, labeled HFB-18, with such a generalized Skyrme force 13]. The 
parameters t$, X5 and 7 were chosen in order to avoid a ferromagnetic collapse of neutron- 
star matter. For simplicity, the remaining parameters in Eq. ( 16 2 p were fixed by the equations 



= 7 , 



X4 



4 + 5x 5 

5 + 4x 5 



(63a) 
(63b) 



(63c) 
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which ensure that the contributions of the new terms to the coefficients CJ vanish identically. 
As a result, the t± and £5 terms cancel exactly in unpolarized homogeneous nuclear matter. 
This new model yields almost as good a mass fit as our previous model HFB-17, with the 
advantage that NeuM matter is now stable with respect to any degree of spin polarizations. 
Even thouth this new force still predicts an isospin instability in SNM, this does not affect 
the interior of neutron stars which is now unpolarized. Moreover we have found that this 
isospin instability can be easily removed if the conditions ( I63aj) — (163c[) are released, without 
deteriorating the quality of the mass fit [50|. However we did not succeed in constructing 
a nuclear mass model that satisfies Eq. ( H7|) . As a consequence, nuclear matter could still 
become unstable at finite temperatures even though no phase transitions occur at zero 



temperature, as shown in Ref. 



One might be tempted to enforce the stability conditions Cj = by adding a zero- 
range tensor force to the conventional Skyrme interaction ([5]) with suitable adjustments of 



the parameters, like the parametrization T22 of Ref. 44j. Unfortunately, a tensor force 



introduces new terms in the functional which also affect the stability of nuclear matter [5 
The stability of 41 different Skyrme interactions having a tensor component has been recently 
studied in Ref. [58]. In particular, the recent Skyrme forces from the Saclay-Lyon group [44J 
which include tensor forces and which were fitted following the same protocol as the older 

n 

SLy family [35J], still predict various spin and spin-isospin instabilities. This is notably the 
case for the force T22 for which Cj = 0. 



V. CONCLUSION 



Nuclear energy density functional theory has been traditionally restricted to very specific 
ogical semi- local functionals of the form given by Eqs. ( |2al) -(l2cl). based on effec- 



Eft 



phenomeno 

tive forces tl 01 . However the use of effective forces introduces tight correlations between 
different terms of the functional, which can generate various kinds of instabilities. In partic- 
ular, the time-odd terms St ■ Tt induced by the momentum-dependent part of Skyrme forces 
(which contribute also to the coupling constants C t T , Cf p and C t As ) are responsible for spu- 
rious spin and spin-isospin instabilities in infinite homogeneous nuclear matter at densities 
encountered in the interior of neutron stars. In some cases, instabilities arise in symmetric 
nuclear matter below saturation densities and could thus also contaminate calculations in 



21 



finite nuclei (Note that the coupling constants C f Ap alone were found to drive finite-size 
59[). These correlations between different parts of the nuclear energy density 



instabilities 



functional hamper the development of more accurate functionals since adding one term in 
the effective force can induce several new terms in the functional. Moreover, the coupling 
constants of the time-odd terms are generally not directly fitted to experimental data but are 
calculated a posteriori using the parameters of Skyrme force. However, there is no guarantee 
that the effects associated with the time-odd terms will be correctly described in this way. As 
shown in Refs. jj, 8| it is not possible to avoid spurious transitions to spin-ordered phases in 
nuclear matter above 2-3 times saturation density. The critical densities above which these 
instabilities occur, decrease when more nuclear data are included in the fit of the parame- 
ters of the SXyrnre force Q. For .stance, for our nuclear mass .node. HFB-17 



40 



41), the 



ground state of neutron matter becomes ferromagnetic above 0.17 fm 3 . These instabilities 



can be (at least partial 



for instance in Refs. 



.y) removed by suitable extensions of the Skyrme force, as proposed 



12 



13| . However an unphysical spin-ordering could still occur at finite 



temperatures thus spoiling the application of Skyrme forces to the hot nuclear matter found 
in protoneutron stars and supernova cores. Alternatively the terms St ■ Tt that are respon- 
sible for spin and spin-isospin instabilities could be canceled by suitable adjustments of an 
additional tensor component to the Skyrme force 44|. Unfortunately a tensor force would 



also generate new terms in the energy density which still lead to instabilities 

On the other hand, the concept of effective forces leads to formal inconsistencies as 
recently discussed in Ref. 60J. Lots of efforts are now devoted to the construction of 
non-empirical functionals from realistic interactions directly without resorting to effective 



forces 



32| . If one adopts the point of view that the nuclear functional is more fundamen- 



tal than effective forces, the different terms appearing in Eqs. (l2b|) and (!2c|) can be treated 
independently (apart from the requirements of gauge invariance and cancellation of self- 
interactions). It is therefore perfectly legitimate to set C/ = — Cj = 0. Actually the J 2 and 
J 2 terms are dropped in most Skyrme forces, not only because of simplicity but also because 



it seems to be favored by global fits to nuclear data and basic nuclear matter properties 611 ]. 
Moreover, the J 2 and J 2 terms might even lead to instabilities in the single-particle spectra 
of finite nuclei, as discussed for instance in Ref. (44 However in all previous studies of spin 
and spin-isospin instabilities in nuclear matter [7 



14,07 



the associated time-odd terms 



s ■ T and (s n — S- 



T p ) have been included in the residual interaction thus violating 
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gauge symmetry. We have therefore reexamined the stability of nuclear matter by setting 
Cf = for those Skyrme parametrizations which were fitted without the J 2 and J 2 terms. 
We have found that this simple prescription not only improves the values of the Landau 
parameters Go, G' Q and Go euM . But this also generally removes all kinds of instabilities in 
asymmetric nuclear matter both at zero and finite temperatures. Nevertheless this prescrip- 
tion yields unrealistic values of the Landau parameters G%, G[ and Gf euM , hence also of the 
effective masses M* a in polarized matter. Further improvements thus require extensions of 
existing Skyrme functionals. 
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TABLE I: Landau parameters Go and G' in symmetric nuclear matter and Gq in neutron matter 
(at saturation density) for selected Skyrme forces which were fitted without the J 2 and J 2 terms. 
Values in parenthesis were obtained by setting Cj = 0. The last line shows the Landau parameters 



predicted by microscopic calculations using realistic interactions: Ref. 
matter and Ref. 1491 for neutron matter. 



481 ] for symmetric nuclear 





Go 


Gq 


Org 


SGII 


0.01 (0.62) 


0.51 (0.93) 


-0.07 (1.19) 


SLy4 


1.11 (1.39) 


-0.13 (0.90) 


0.11 (1.27) 


Skll 


-8.74 (1.09) 


3.17 (0.90) 


-5.57 (1.10) 


SkI2 


-1.18 (1.35) 


0.77 (0.90) 


-1.08 (1.24) 


SkI3 


0.57 (1.90) 


0.20 (0.85) 


-0.19 (1.35) 


SkI4 


-2.81 (1.77) 


1.38 (0.88) 


-2.03 (1.40) 


SkI5 


0.28 (1.79) 


0.30 (0.85) 


-0.31 (1.30) 


SkO 


-4.08 (0.48) 


1.61 (0.98) 


-3.17 (0.97) 


LNS 


0.83 (0.32) 


0.14 (0.92) 


0.59 (0.91) 


Realistic 


0.83 


1.22 


0.77 
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TABLE II: Critical densities above which nuclear matter becomes unstable according to Landau's 
criterion for selected Skyrme forces which were fitted without the J 2 and J 2 terms. The first two 
column are for symmetric nuclear matter, while the last column is for pure neutron matter. The 
densities indicated in parenthesis were obtained by setting Cf = 0. 





p c (G ) [fm- 3 ] 


Pc(G' Q ) [fm- 3 ] 


Pc (G^ M ) [fm- 3 ] 


SGII 


0.44 (oo) 


0.80 (oo) 


0.26 (2.07) 


SLy4 


oo (oo) 


0.33 (oo) 


0.59 (oo) 


Skll 


0.04 (0.71) 


oo (oo) 


0.05 (oo) 


SkI2 


0.14 (oo) 


oo (oo) 


0.15 (oo) 


SkI3 


0.91 (oo) 


0.92 (oo) 


0.37 (oo) 


SkI4 


0.07 (oo) 


oo (oo) 


0.09 (oo) 


SkI5 


0.43 (oo) 


1.36 (oo) 


0.28 (oo) 


SkO 


0.07 (0.52) 


oo (2.32) 


0.09 (0.67) 


LNS 


oo (oo) 


0.43 (oo) 


0.62 (1.38) 



TABLE III: Landau sum rules given by Eqs. (|59ap and (|59b|) for selected Skyrme forces which were 
fitted without the J 2 and J 2 terms. Values in parenthesis were obtained by setting Cj = 0. 





Si 


s 2 


SGII 


0.97 (0.61) 


1.13 (-0.51) 


SLy4 


-0.31 (-0.65) 


1.52 (0.85) 


Skll 


-6.71 (-0.59) 


-89.2 (0.86) 


SkI2 


6.87 (-0.71) 


-20.7 (0.98) 


SkI3 


-1.46 (-2.33) 


2.14 (1.84) 


SkI4 


1.01 (-1.23) 


-11.3 (1.32) 


SkI5 


-1.47 (-2.28) 


2.17 (1.77) 


SkO 


3.21 (1.07) 


-13.7 (0.87) 


LNS 


0.49 (0.63) 


3.53 (-0.04) 
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TABLE IV: Landau parameters Go an d G' in symmetric nuclear matter and Gq e in neutron 
matter (at saturation density) for selected Skyrme forces which were fitted with the J 2 and J 2 
terms. The last line shows the Landau parameters predicted by microscopic calculations using 
realistic interactions: Ref. 



48j for symmetric nuclear matter and Ref. 



49j for neutron matter. 





Go 


G 


/^yNcuM 
G 


SkO' 


-1.62 


0.79 


-1.43 


SLy5 


1.09 


-0.16 


0.09 


SkP 


-0.23 


0.06 


-0.61 


SkX 


-0.63 


0.51 


-0.50 


BSkl7 


-0.69 


0.50 


-0.88 


BSkl7st 


-0.68 


0.50 


0.47 


BSkl8 


-0.33 


0.46 


-0.57 


Realistic 


0.83 


1.22 


0.77 



TABLE V: Critical densities above which nuclear matter becomes unstable according to Landau's 
criterion for selected Skyrme forces which were fitted with the J 2 and J 2 terms. The first two 
column are for symmetric nuclear matter, while the last column is for pure neutron matter. 





p c (G ) [fm- 3 ] 


p c (G' ) [fm" 3 ] 


Pc(G^ M ) [fm- 3 ] 


SkO' 


0.12 


0.97 


0.14 


SLy5 


oo 


0.33 


0.57 


SkP 


0.74 


0.30 


0.19 


SkX 


0.22 


0.40 


0.19 


BSkl7 


0.21 


0.68 


0.17 


BSkl7st 


oo 


oo 


oo 


BSkl8 


oo 


0.62 


oo 
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